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We study the geometric phase factors underlying the classical and the corresponding quantum 
dynamics of a driven nonlinear oscillator exhibiting chaotic dynamics. For the classical problem, we 
compute the geometric phase factors associated with the phase space trajectories using Frenet-Serret 
formulation. For the corresponding quantum problem, the geometric phase associated with the time 
evolution of the wave function is computed. Our studies suggest that the classical geometric phase 
may be related to the the difference in the quantum geometric phases between two neighboring 
eigenstates. 
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Since the discovery of Berry phase T|, the study of geo- 
metric phases based on the common mathematical theme 
of anholonomy, has emerged in a variety of fields. 0] The 
Berry phase is a path-dependent geometric phase associ- 
ated with the adiabatic time evolution of the wave func- 
tion, associated with circuits iiiparameter space. This 
concept has been extendedy,!!,!^] to non-adiabatic cases 
and also to non-cyclic circuits, since the phase acquired 
by the wave function in any type of time evolution, may 
have a component that is of purely geometric origin. This 
phase is a gauge invariant quantity and is equal to the dif- 
ference between the total phase and the dynamical phase 
acquired by the wave function. 

One of the open questions has been the classical limit 
of the Berry phase or its generalization describing the 
geometric part of the phase of the quantum wave func- 
tion. For the special case of integrable Hamiltonian sys- 
tems, described in terms of action-angle variables, the 
so-called Hannay angle BfJ^ represents the semiclassical 
limit of the Berry phase. Berry gave an explicit formula 
relating the classical angle, called the Hannay angle and 



the n quantum eigenstate (/>„ as, 6^ 
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There 



have been some attempts to describe the classical limit 
of Berry phase for chaotic systems where the effort has 
been focused on finding a generalization of the 2-form 
within Wigner-Weyl formalism. Q 

Viewing the Berry phase as an anholonomy effect un- 
derlying dynamical evolution described by Schrodinger 
equation, we seek a classical analog of anholonomy, un- 
derlying the corresponding classical evolution described 
by the Newton's equations of motion. Here, we are con- 
cerned with the geometric phases associated with pe- 
riodic, quasiperiodic and chaotic dynamics of a driven 
nonlinear oscillator. We compute the geometric com- 
ponent of the phase of the wave function using a kine- 
matic formulation^ of the Berry phase Q| as given by 
Mukunda and Simon. In the corresponding classical 
problem, we find the anholonomy underlying nonplanar 



periodic phase space trajectories and then extend this 
formulation to quasiperiodic and chaotic trajectories. It 
should be emphasized that unlike in Berry phase, our cir- 
cuits are not in parameter space but are in phase space. 
By treating a classical phase space trajectory as a space 
curve, we show that a geometric phase can be associ- 
ated with every trajectory, by using a Frenet-Serret (FS) 
formulation 8|. The classical geometric phase is the in- 
tegrated torsion of the phase trajectory, and it bears a 
strong analogy to the geometric phase factor associated 
with the wave function. 

To illustrate this correspondence, we begin with the 
Frenet-Serret equations, describing the time evolution of 
the orthonormal FS triad, consisting of the tangent T, 
the normal N, and the binormal B to the space curve 



vkN . N : 



-vkT + vtB,B 



-vtN . 
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where k and r are respectively the curvature and tor- 
sion of the curve and v = |r|. Here, the overdot denotes 
derivative with respect to time. The above equations can 
be written as, F = ^ x F, where F = T, N, or B, and 
^ = —V kB -I- vtT. That is, the FS triad rotates around 
T and B. One way to quantify this rotation is to work in 
a frame in which T is parallel transported, and measure 
the angle of rotation around the tangent T. This will be 
given by the angle (j)c{t) — /q rvdt' . Thus (j)c{t) is the ge- 
ometric phase characterizing the anholonomy associated 
with the corresponding phase space orbit. 

If we define a complex vector M = (N -I- iB) /V2, the 
classical geometric phase can be written as 
(j>^{t) = i/J M* • M vdt'. 

This expression has exactly the same form as the quan- 
tum geometric phase found by Berry [J , when the clas- 
sical vector M is replaced by a quantum eigenstate ipn- 
Additionally, by mapping the closed phase space trajec- 
tory to a circuit on a unit sphere traced by the tip of the 
tangent vector ( tangent indicatrix), (j)c can be shown to 
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be the solid angle subtende d by this tangent indicatrix 
at the center of the sphere. These results are valid 
also for a non-periodic trajectory, since it can always be 
closed using a geodesic on the sphere. 4] 

Motivated by this close analogy between the geometric 
phase in a quantum system and the FS geometric phase, 
we investigate any possible relationship between the two. 
The underlying key question is whether the classical an- 
holonomy is related to the corresponding quantum one. 
Here we calculate the classical and the quantum geomet- 
ric phases for a periodically driven nonlinear oscillator 
exhibiting complex dynamics. The system under inves- 
tigation is an "impact oscillator"!^, the oscillator that 
rebounds elastically whenever its displacement x drops 
to zero. For a; > 0, the system is described by, 

H=^p^ + ^Lj^x^ - / cosiu;t)x (2) 

The system is piecewise linear, and the analytic solutions 
can be obtained for a; > 0. The discontinuity at the origin 
makes it essentially nonlinear. Without loss of generality, 
we choose the units of t and x such that uj = 1 and 
/ = 1. The phase space trajectory of the dynamical 
system can be viewed as a space curve generated by the 
three-dimensional vector r(t) ~ (x,x,x) parameterized 
by the time t. 

As we discuss below, the impact oscillator exhibits very 
rich and complex dynamics, where periodic, quasiperi- 
odic and chaotic dynamics coexist. Our choice of this 
example was motivated by the fact that in addition to 
the simplicity underlying the classical analysis of the os- 
cillator, the quantum wave functions for the driven os- 
cillator are known in terms of the classical solution as, 

i>ix, t) = xix', t) exp + L{t')dt'] (3) 

Here x' = x — Xc{t), with Xc{t) being the solution of 
the classical equation of motion and L is the Lagrangian 
of the driven system, xi^i t) is the wave function of the 
oscillator in the absence of driving. Note that the wave 
functions of the driven oscillator are centered on the po- 
sition of the classical forced oscillator. 

If we take x to be an eigenstate of the undriven os- 
cillator, Xn{x,t) = Un{x) exjp —i{Ent / h) , the eigenfunc- 
tion can be written in terms of Hermite polynomials as 
Un{x) — exp~[u!o/{2h)x'^]Hn{\/uJo/hx). It should be 
noted that corresponding to the eigenstate, \'tp{x,t)\'^ — 
\Xn{x' jt)]"^ . Therefore, the center of the wave packet 
Xc{t) obeys classical equation of motion, and the shape 
of the wave packet ( the density distribution with respect 
to the center Xc) is unaffected by the driving force. 

In view of the fact that we have the analytic solution 
for the classical system for x > 0, and also a closed form 
solution for the quantum wave function, we can compute 



the classical and the quantum geometric phases with ex- 
treme precision. These will be presented below. 

Figures 1 shows richness and complexity underlying 
the classical dynamics of the oscillator as we vary the 
initial energy (initial conditions) of the oscillator. We 
see periodic orbits and quasipcriodic tori, in addition to 
chaotic trajectories. 

All results are for a fixed loq = 1.6 which corresponds to 
the oscillator frequency of 3.2, in units of w, the frequency 
of the driver, putting our analysis is close to the adiabatic 
regime. 

Furthermore, we believe that the probability of induc- 
ing a transition (due to driving) is rather small because 
the classical energies of the particle under consideration 
here are far below the threshold for the transition be- 
tween two quantum states ni and n2 in units of Tlloq. 
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FIG. 1: For po = 0, the figure shows the possible x and p 
values (once every period of the driver), as a function of xq, 
the initial position of the oscillator. For xo ~ 1.0416, there is 
a period- 11 orbit. For xq > 1.046, we get invariant quasiperi- 
odic tori or cantori as seen in the lower figure. For xq < 1.035, 
almost all initial conditions result in chaotic dynamics. 

Using the classical solution, the FS geometric phase 
(j)c can be computed as the integral of the torsion, given 
by r = rt • {ra x rttt)/\rt x rttp For the corresponding 
quantum problem, we can obtain the geometric phase, 
using kinematic formulation by Mukunda and SimonQ. 
For a given wave function i/;(a;, i), the quantum geometric 
phase 4'qis a gauge invariant quantity is given by (j)q{t) = 
arg J[ip*{x,0)ip{x,t)]dx - Im J^^ ip* {x,t')dtil^{x,t')dt' . It 
is clear that the first term describes total phase while the 
second term describes the dynamical phase accumulated 
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by the wave function in time t. Note that this formulation 
does not require any circuits to define geometric phase. 

For the driven impact oscillator, if we substitute the 
solution for the wave function, given by equation (3). , 
the geometric phase can be written in terms of the ex- 
pectation values of < a; > and the classical solution Xc 



- n<j>q.n = / [L{t')- <x> Xcit')]dt' + G{t) (4) 
Jo 

where G{t) is given by, G{t) = (xc{t)xc{t) — Xc{Q)xc{Q)) + 
arg[j[{unix'{0)unix'{t) exp 'f{xc{t) - Xc{0))]dx. For pe- 
riodic evolution, G{t) — 0. Also, G{t) = if we consider t 
values where the classical particle is at the turning points. 
In view of this, we will confine our calculations of geo- 
metric phases to only such values of t. 

Since the classical impact oscillator is confined to a; > 
0, the eigenstates of the corresponding quantum system 
are restricted to odd n values only. Substituting the ex- 
plicit expression for the Hermite polynomials, < x > can 
be expressed in terms of parabolic cylindrical functions 
which are functions of Xr- 



< X >2«-1 = 



(5) 



where D-miUc) is a parabolic cylinder functions of order 
771. Here i/c — axc where a = ^/2uJoJh. Ai{yc) and 
Bi{yc) are polynomials of yc of degree I that are uniquely 
determined by n. For example, for n = 1, we have 



1 6i?_4(yc) - ^ycD-3{y,) + ylD^2iyc) 



iiyc) - 2y,D.2{yc) + ylD_^{y,) 



(6) 



For higher values of n, the expressions are very com- 
plicated. Our analysis has been confined to the ground 
state n = 1, and the first excited state n = 3. For all val- 
ues of n, Ti4>q,n reduces to the classical action as ?i ^ 0. 
We factor out this part of the phase factor 0o and write 
4>q,n = 4>o + 'Pn- As we show below, it is the difference in 
the phases between the two neighboring eigenstates that 
exhibit quantum fingerprints of classical dynamics. 

Figures 2 and 3 show geometric phases for a fixed time 
interval T , equal to the driving period, for a periodic and 
a chaotic trajectory. These results as well as other similar 
analysis suggest that (pc may be related to (f>n~i — 0n- 

For dynamical evolution involving arbitrary time t, 
our detailed analysis shows that the integrated quantum 
phase factors can be expressed as a sum of linear and 
oscillatory functions of t and hence can be written as , 
<pn(t) = Vnt + fn{t). Hcrc v„ are constants, independent 
of t and fn{t) are oscillatory functions which are found 
to retain the fingerprints of the corresponding classical 
dynamics for all times. 

In contrast to quantum phases, the classical 4>c{t) 
is found to be an oscillatory function for periodic , 



quasiperiodic as well as for the chaotic trajectories. This 
should be contrasted with the the driven and damped 
oscillator phases 13j where the classical phase averaged 
over the period of the driver was finite. 




FIG. 2: For a period-11 trajectory, the four figures (from 
top to bottom) show the quantum geometric phases for the 
ground state( n — 1), the first excited state { n — 3), classical 
geometric phase (/f>c and dcj> = cj>i — 4>3, respectively. 

The results for arbirary time evolution are shown in 
figures 4 and 5. Here we factor out the linearly increas- 
ing parts of the phases and show the fluctuating parts 
along with the classical phase. For the parameter values 
corresponding to the figures 2 to 5, vi/{2tt) — .4218 and 
W3/(27r) = .3956 for the periodic orbit and vi/{2-k) — 
.3723 and W3/(27r) = .4021 for the chaotic orbits. We no- 
tice that 4>c is of the same order of magnitude as /i — /a 
with (pc exhibiting intermittent fluctuations. In view of 
the fact that vi ~ fs, it is conceivable that Vn approaches 
a constant, independent of n as n ^ oo and therefore (pc 
may be related to the for arbitrary time t. This is 
analogous to Berry's relation 0„ = —dntpn which was 
true for integrable systems in semiclassical limit. 

In summary, our results as depicted above describe pre- 
liminary studies of classical and quantum phases under- 
lying a driven nonlinear oscillator. An interesting result 
is the possible relationship between the classical phases 
and the difference in the quantum phases between the 
two neighboring states, reminiscent of the Berry relation 
6h — —dn4>n- As a consequence, the smalless of the clas- 
sical phase will have its origin in the relative stability of 
the quantum phase with respect to the quantum num- 
ber n. It is rather intriguing that the fluctuations in the 
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FIG. 3: For a chaotic trajectory, same plots as in Fig. 2 In 
the third figure, we superimpose dcf) with 0c. 
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FIG. 4: For a poriod-11 trajectory, the top two figures show 
the oscillatory components /i and /a of the net accumulated 
phases for the ground state and the excited state. The third 
figure shows the classical phase and the differnce d/ = /i — /a . 
The lowest figure shows the position of the oscillator. 



FIG. 5: For a chaotic trajectory, the same plots as in Fig. 4. 



quantum geometric phases retain the finger prints of the 
corresponding classical dynamics. Further detailed stud- 
ies involving higher excited states are needed to confirm 
these speculative views. 

In physical applications such as atom optics, Hamil- 
tonians of the form, H{x,t) = Hq + Xx sm(ojt) are rele- 
vant where Hq is the time-independent Hamiltonian and 
the time-dependent term describes the interaction with 
a single mode radiation field in dipole approximation. 
For nonlinear Hq, nontrivial dynamics may lead to many 
surprises. The impact oscillators exhibit dynamics with 
many features that are typical of a nonlinear oscillator. 
This suggests that the driven impact oscillator may pro- 
vide an interesting theoretical model to explore various 
issues relevant to quantum chaos. 
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